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^ The attractive branch of the interaction potentials with the surface of gold 

. ^ have been computed for a large variety of atomic systems: the hydrogen atom, 

noble gases (He, Ne, Ar, Kr, Xe), alkali atoms (Li, Na, K, Rb, Cs) and alka- 
line atoms (Be, Mg, Ca, Sr, Ba). The results include highly accurate dynamic 
polarizabilities for the helium atom calculated using a variational method and 

O-h explicitly correlated wavefunctions. For other atoms considered we used the 

data available in the literature. The interaction potentials include both the 
effects of retardation of the electromagnetic interactions and a realistic repre- 

^ sentation of the optical response function of gold (beyond the approximation 

"Zl of a perfect conductor). An explicit comparison of our result to the interaction 

between an atom and a perfect conductor is given. 
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1. Motivation 

m 

t-H In the last decade, there has been drastic experimental progress in the field 

of (ultra-) cold collisions in atomic systems. These may involve interac- 

J^j tions among two atoms as well as those between an atom and a solid wall. 

Among the latter, one of the experimental methods for studying the long 
range (z ^> cto), attractive part of the interaction potential makes use of 
the phenomenon of quantum reflection. It was demonstrated, that scatter- 
ing cold atomic beams under grazing angles from well defined, single crystal 
solid surfaces is especially useful for probing Casimir-Polder forces. This is 
due to the fact that quantum reflectivity is particularly sensitive to the 
shape of the potential, at distances where it is heavily modified by retarda- 
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tion. 1 Other techniques involving the manipulation of clouds of cold atomic 
gases, or condensates in the presence of a solid wall 2 heavily depend on the 
details of the interaction potential. Despite the abundance of experimental 
results, the theoretical analysis of the interaction potentials was often based 
on simple model potentials with parameters chosen to fit the data. Here, 
we provide accurate atom-surface interaction potentials for the surface of 
gold based on ab-initio computed atomic dynamic polarizabilitics and a 
compilation of all available optical data in literature. 

2. Atom-surface interactions 

For distances where the exchange effects become negligible, the atom- 
surface can be computed by considering a polarizable particle interacting 
with quantum electromagnetic field fluctuations. In this approach, the solid 
is treated as a continuous medium having a frequency dependent permit- 
tivity e{uj). The derivation of the dipolar term V\(z), dominant at the long- 
distance limit, has been first performed by Lifshitz, 3 and the result when 
given in atomic units (ra e = e = h = 1) reads: 

3 poo />oo 

Vi{z) = -%- oWaiM/ d£c- 2 ^H(£,6(ic)) , (i) 
27r Jo J i 

where a is the fine structure constant and ai(iw) is the dipole polarizability 

of the atom [for the fundamental physical constants we take the COD ATA 4 

recommended values, e.g., a = 137.035999679(94)]. The function H(£,e) is 

given by: 

m e) = a - 2a + > f+^ -* . (2) 

This expression simplifies considerably in the limit of a perfect conductor 
[e(uj) — > oo], for which case the potential becomes: 

V}°° ) (z) = / &Lu ai (iLo)c- 2a " z P^iauz) , (3) 

where p[ ca \z) = I + 2z + 2z 2 . Another simplification takes place in the 
short-distance and long-distance limits. For z — > 0, the interaction potential 
behaves as: 

VAz) — -- — 7T / dwaWiw) — — - = 5-. (4) 

The long-distance limit the interaction potential for a generic e(uj) is: 

3ai(0) 6(0) -1 C\ 
Vl(z) ~ ~ W e(0) + l=-^- (5) 
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but this result cannot be used for the case of conductors, where e has a 
pole at oo=0. We have checked that when the e(iw) (a real function) is 
bounded from below by a/ui, as it is in the case of conductors, the long 
distance behavior of the potential is equal to the one of perfect conductor, 
asymptotically parametrized by C4 00 ': 

Vr nd (z) z ~°° ^ (oo) (z) z - 3ai( -°) = - (6) 
1 w ivy 87raz4 z4 v J 



3. Frequency dependent dielectric permittivity of gold 

The ab-initio computation of the frequency dependent dielectric response 
function of metals is beyond the reach of present day electronic structure 
calculations. In this work, the complex frequency dependent permittivity of 
gold was reconstructed from the experimental optical data. The imaginary 
part of the dielectric constant was modeled using a function: 

Ime(w) = p T +Vc n f T LK„,wi n ,7n;w) , (7) 
lj(oj' ! +oj£) ' 

v T 1 n 

where the first term is the Drude model of the free-electron contribution, 
and the second is a sum of empirical Tauc-Lorentz 5 functions: 

where 6{x) is the Heaviside step function [8(x) = for x < 0, and 6(x) = 
otherwise]. Once a satisfactory representation of the imaginary part of 
e(uj) is found, the real part of the permittivity can be calculated using the 
Kramers-Kroenig relation. The parameters of Eq. (7) have been fitted to the 
experimental optical data. The data used included the data sets collected 
in the handbook by Palik 6 which cover the visible, ultraviolet and X-ray 
regimes, and various compilations of optical data in the microwave and in 
the infrared regions of the electromagnetic spectrum. 7 

The results of a global fit of the model (7) to the experimental data for 
gold are listed (in atomic units) in Table 1 and depicted in Fig 1. Our fitted 
values of the plasma frequency (w p ), and the damping frequency (oj t ) can 
be compared to the values used by Lambrecht et al. s who used values of 
w p =0.330, w T =0.00108 and cj p =0.276, w T =0.00478. 

When using the model defined in Eq. (7), we observe a slight, but sys- 
tematic deviation of the fitted function with respect to the experimental 
data for the imaginary part of the permittivity in the very low frequency 
region. In order to improve the fit, we decided to further refine our model 
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Table 1. Best fit parameters for the analytic model of Ime(u) for 
gold according to Eq. (7). 
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0.001636 
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"5 
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75 
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for the dielectric function by including an additional Drude-model term, 
corresponding to a conductor with two types of carriers: 



,2 , , , ,2 



The parameter of the model function have been reoptimized using nonlin- 
ear least squares fit to the optical data, and their values are presented in 
Table 2. The double Drudc model leads to a significant improvement of the 



Table 2. Best fit parameters for the two plasma frequency model of 
Imt(u) for gold according to Eq. (9). 
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low-frequency behavior of e, but has a negligible influence on the calculated 
potential. Since there is no physical justification for Eq. (9), we use Eq. (7) 
for the results obtained in the rest of this paper. 



4. Numerical Results 

The atom-surface interaction potentials have been computed by numeri- 
cally evaluating Eq. (1) using Eqs. (7) and (9) and ab-initio calculated 
values of dynamic dipolc polarizabilities. For the case of helium the re- 
spected polarizabilities have been computed variationally using basis sets 
of explicitly correlated functions. 9 For the other atoms considered we used 
dynamic polarizabilities published by Derevianko, Porsev and Babb. 10 
We find it convenient to represent the potential by its short distance 
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limit multiplied by a "damping function" , accounting for the effects of re- 
tardation, which for the dipole case reads: 

Vi(z) = -^Mz), (10) 

For the damping function fs(z), we found the following functional form to 
lead to an accurate representation of all interaction potential for all atom- 
wall distances: 

f ( \_ 1 + «i az + a 2 {az) 2 

MZ > ~ l + ai az + b 2 (az) 2 + b 3 (az) 3 ' [ ' 

with a being the fine structure constant. The equality of the linear terms 
in both the numerator and in the denominator of this rational function is 
due to the requirement of the potential having the correct short distance 
limit: fa{z) = 1 + 0[(az) 2 ]. The best values of the parameters a\, a 2 , b 2 , 
63 in (11) fitted to the calculated interaction potentials for different atomic 
species are presented in Tables 3-5, together with the corresponding values 
of C 3 . 



Table 3. Long range potentials for noble gas atoms interacting 
with a surface of gold. The second row contains the C3 constants, 
and the following rows give the values of the /3(z) damping func- 
tion, defined in Eq. (10) for different distances. In the last four 
rows, we present the best fit parameters for the rational func- 
tion (11) which covers all distances. 





He 


Ne 


Ar 


Kr 


Xc 


Cz 


0.062(3) 


0.127(1) 


0.415(4) 


0.588(5) 


0.870(6) 


z 




AW 


= V(z)/(C 3 / 


z A ) 




lxlO 1 


0.99108 


0.98804 


0.99323 


0.99348 


0.99405 


2X10 1 


0.97388 


0.96653 


0.98014 


0.98112 


0.98275 


5X10 1 


0.90949 


0.89155 


0.92931 


0.93352 


0.93907 


lxlO 2 


0.80490 


0.77796 


0.84192 


0.85135 


0.86320 


2xl0 2 


0.64602 


0.61531 


0.69933 


0.71522 


0.73546 


5xl0 2 


0.40337 


0.37942 


0.45908 


0.47873 


0.50522 


lxlO 3 


0.25046 


0.23481 


0.29294 


0.30931 


0.33232 


2xl0 3 


0.14492 


0.13575 


0.17185 


0.18274 


0.19849 


5xl0 3 


0.06539 


0.06124 


0.07800 


0.08322 


0.09088 


lxlO 4 


0.03441 


0.03222 


0.04109 


0.04387 


0.04796 


2xl0 4 


0.01770 


0.01657 


0.02114 


0.02257 


0.02469 


5xl0 4 


0.00721 


0.00675 


0.00861 


0.00919 


0.01005 


lxlO 5 


0.00363 


0.00339 


0.00433 


0.00462 


0.00506 


ai 


4.19366 


5.91822 


3.76881 


3.95140 


4.21331 


(12 


0.29208 


0.47236 


0.21430 


0.20804 


0.20365 


b 2 


2.13922 


3.33497 


1.53023 


1.47551 


1.41033 


b 3 


0.10861 


0.18758 


0.06671 


0.06066 


0.05430 
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Tabic 4. Long range potentials for alkali atoms interacting with 
a surface of gold. The second row contains the C'z constants, and 
the following rows give the values of the fz(z) damping function, 
defined in Eq. (10) for different distances. In the last four rows, we 
present the best fit parameters for the rational function (11) which 
covers all distances. 





Li 


Na 


K 


Rb 


Cs 


Cz 


1.210(5) 


1.356(6) 


2.058(9) 


2.31(1) 


2.79(1) 


z 






= V(z)/(Cz/ 






lxlO 1 


0.99899 


0.99845 


0.99829 


0.99801 


0.99780 


2X10 1 


0.99740 


0.99590 


0.99542 


0.99465 


0.99408 


5X10 1 
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0.98533 


0.98333 


0.98086 
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lxlO 2 


0.97415 


0.96515 
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0.95541 
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2xl0 2 


0.94032 


0.92509 


0.91803 


0.90850 


0.90210 


5xl0 2 


0.84866 


0.82359 


0.82012 


0.80456 


0.79521 


lxlO 3 


0.73176 


0.69946 


0.70857 


0.69126 


0.68309 


2xl0 3 


0.57369 


0.53721 


0.56295 


0.54776 


0.54485 


5xl0 3 


0.34082 


0.30970 


0.34388 


0.33492 


0.33918 


lxlO 4 


0.19776 


0.17688 


0.20301 


0.19808 


0.20323 


2xl0 4 


0.10552 


0.09370 


0.10926 


0.10672 


0.11029 


5xl0 4 


0.04353 


0.03854 


0.04525 


0.04422 


0.04585 


lxlO 5 


0.02196 


0.01943 


0.02285 


0.02233 


0.02318 


Ol 


6.96178 


19.81680 


1.50957 


1.90225 


2.21418 


0-2 


0.10676 


0.30362 


2.97205 


2.06782 


1.55523 


*>2 


0.44672 


1.43152 


3.15260 


2.30121 


1.83158 


bz 


0.00659 


0.02118 


0.15449 


0.11044 


0.07999 



The potentials presented above can be compared to simplified ones used 
previously 1 in the analysis of the atomic beam experiment measuring the 
quantum reflectivity There, retardation has been accounted for using: 

V < S) = -?(7TA)' (12 » 

in which A represents the reduced wavelength of the first electronic dipole 
transition in He, i.e. A = 178 Bohr (93A). This function has also been used in 
the preliminary evaluation of the experimental results for 3 He atoms quan- 
tum reflecting from a single crystal gold surface. Here however, significant 
deviations have been found indicating, that the coarse-grained model (12) 
is not sufficiently complex to mimic the more involved rational structure 
given in Eqs. (10) and (11). 

An independent check of our potentials has been performed by using 
our model for the dielectic constant e(iw) of gold to calculate the interac- 
tion energy between two surfaces of gold, for both semi-infinite solids and 
for thin foils. Our computations are within 1% of the revisited results of 
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Table 5. Long range potentials for alkaline atoms interacting with 
a surface of gold. The second row contains the C3 constants, and 
the following rows give the values of the fs(z) damping function, 
defined in Eq. (10) for different distances. In the last four rows, we 
present the best fit parameters for the rational function (11) which 
covers all distances. 





Be 


Mg 


Ca 


Sr 


Ba 




0.650(5) 


1.067(6) 


1.82(1) 


2.19(1) 


2.75(1) 


z 




h(z) 


= V(z)/(C 3 / 






lxlO 1 


0.99797 


0.99790 


0.99805 


0.99788 


0.99778 


2X10 1 


0.99388 


0.99394 


0.99442 


0.99404 


0.99382 


5X10 1 


0.97571 


0.97722 


0.97909 


0.97810 


0.97764 


lxlO 2 


0.93887 


0.94441 


0.94957 


0.94792 


0.94755 


2xl0 2 


0.86495 


0.87903 


0.89248 


0.89054 


0.89129 


5xl0 2 


0.69195 


0.72244 


0.75855 


0.75885 


0.76468 


lxlO 3 


0.51642 


0.55531 


0.61130 


0.61588 


0.62880 


2xl0 3 


0.34021 


0.37699 


0.44043 


0.44893 


0.46827 


5xl0 3 


0.16580 


0.18850 


0.23484 


0.24315 


0.26173 


lxlO 4 


0.08878 


0.10168 


0.12963 


0.13511 


0.14765 


2xl0 4 


0.04590 


0.05270 


0.06773 


0.07078 


0.07784 


5xl0 4 


0.01872 


0.02151 


0.02772 


0.02899 


0.03196 


lxlO 5 


0.00942 


0.01082 


0.01395 


0.01460 


0.01610 


01 


3.00255 


3.88303 


12.61180 


3.09874 


2.83515 


0-2 


0.12664 


0.12802 


0.19042 


4.88813 


2.70340 


b 2 


0.51498 


0.55406 


1.27280 


5.20283 


3.03090 


63 


0.01818 


0.01601 


0.01848 


0.41434 


0.20650 



Lambrccht et at 11 ' 12 This gives us confidence, that the ready-to-use inter- 
action potentials for the various atomic species presented in this paper are 
accurate. In addition, it demonstrates that our simple analytic model of the 
permittivity of gold may be very useful for future computation of Casimir 
potential involving gold, for both microscopic and macroscopic bodies. 
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